Fermi-Edge Resonance and Tunneling in Nonequilibrium Electron Gas 
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Fermi-edge singularity changes in a dramatic way in a nonequilibrium system, acquiring features 
which reflect the structure of energy distribution. In particular, it splits into several components if 
the energy distribution exhibits multiple steps. While conventional approaches, such as bosonization, 
fail to describe the nonequilibrium problem, an exact solution for a generic energy distribution can 
be obtained with the help of the method of functional determinants. In the case of a split Fermi 
distribution, while the 'open loop' contribution to Green's function has power law singularities, the 
tunneling density of states profile exhibits broadened peaks centered at Fermi sub-levels. 



Properties of quantum systems can change drastically 
when they are driven out of equilibrium. This is espe- 
cially true for transport in nanodevices [1] , such as quan- 
tum dots and quantum wires, where energy relaxation 
takes place outside the device. Transport in interact- 
ing systems is often difficult to describe by the methods 
developed for analyzing equilibrium [2-7], which makes 
exact solutions outside equilibrium scarce and valuable. 

Fermi-edge singularity [8,9] (FES) is a dramatic man- 
ifestation of interactions and correlations in electron liq- 
uid. It can be observed in a resonant tunneling experi- 
ment [10,11] as a power law resonance which peaks at the 
Fermi level. Being one of the few exactly solvable prob- 
lems describing transport in strongly interacting systems, 
FES has been thoroughly explored in a variety of situ- 
ations, including quantum wires [12-14], quantum Hall 
edge states [15] and quantum dots [16]. However, apart 
from recent work [17] which resolved a long-standing con- 
troversy on orthogonality catastrophe in two Fermi seas 
[18-20], little is known about FES out of equilibrium. 

Nonequilibrium electron states with structured en- 
ergy distribution were demonstrated recently [21] using 
diffusion-cooled nanoscale wires. In a wire short enough 
to allow electrons diffuse out without energy relaxation 
a distribution consisting of two Fermi steps. 



n{e) = (1 — x)nF{e — Mi) + xnpie. - ^2) 



(1) 



with /zi_2 potentials in the leads, was created, imaged us- 
ing tunneling spectroscopy and employed to study energy 
relaxation. Similar approach [22] was used to observe 
splitting of a Kondo resonance in a quantum dot with a 
mixture of two Fermi steps injected in one of the leads. 

Here we study how the FES tunneling density of states 
is modified by nonequilibrium electron distribution, and 
find that it can acquire a rich and complex structure. 
Since FES peaks at the Fermi level, one expects a mul- 
tiple FES peak profile for a multi-step distribution of 
Refs. [21,22], with each FES peak centered around a cor- 
responding Fermi sub-level. While the standard methods 
used to describe FES in equilibrium fail, an exact solu- 
tion can be obtained with the help of a method proposed 
below which allows to extend the FES theory to generic 
nonequilibrium systems. 

Although a variety of methods is available to treat the 



FES problem, applying them outside equilibrium is often 
problematic. The original approach [9], based on resum- 
mation of diagrammatic series, is cumbersome and proves 
difficult to generalize. Thus alternative techniques, most 
notably bosonization [23], have been developed. The 
bosonization approach, however, relies too strongly on 
the assumption of thermodynamic equilibrium, and thus 
cannot be used in our problem. 
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FIG. 1. Fermi-edge resonance splitting (2) for the two-step 
Fermi distribution (1), with the scattering phase S = 0.3. 

The method used in this article is free of such limita- 
tion. It allows to represent the Green's function of an 
electron tunneling in a many-body system in terms of an 
appropriate functional determinant and related quanti- 
ties which are defined in a one-particle Hilbert space. The 
determinant structure accounts in an exact way for all the 
effects of Fermi statistics, as well as for the interaction in 
the final state underpinning the FES phenomenon. Here 
we employ a generalization of the method of Ref. [17] 
recently used in nonseparable mesoscopic FES problem 
[24] , which allows to handle an arbitrary energy distribu- 
tion. After developing general formalism we focus on the 
two-step case (1) and obtain a split FES profile ImG(e) 
in terms of the scattering phase shift 6 (Fig.l), where 



G(e) cx 



l-n(6') 



(e'-Mi)"i(e'-M2)"= 



X D{e - e')de' (2) 



with complex exponents ai — 2(6 — S)/tt, — ^bjix and 



(5 = — In (1 - X e^'^x) 
2i 



(3) 
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The two factors in Eq.(2) correspond, as we will sec, to 
the well known separation [9] of equilibrium FES into the 
'open line' and 'closed loop' contributions. The closed 
loop factor D{e) equals e'^'^ in equilibrium. We 

evaluate D(e) below and find that it describes broadening 
of nonequiHbrium FES, with 7 ~ x{l — x)\ni — H2\S^/''^, 
which can be attributed to a finite effective temperature 
~ J n{l — n)de. The relation ai + a2 = 26 /-k ensures 
agreement with the equilibrium FES exponent. 

Turning to the analysis, the FES Hamiltonian de- 
scribes band electrons interacting with a localized state: 

H = Ho6+S + Hi(l-S+6), Ho,i =^/iya+ap, (4) 

pp' 

where h, 6+ describe the localized state occupation and 
/ipp/^ = CpSppi + V^^^'p? are the single particle operators 
of band electrons scattering on the charged/imchargcd 
state potential V'^^'^\r). Tunneling from the localized 
state is described by the Green's function 



G(r) = tr 



+ (-T)&(-r)p) , 



T > 0, 



(5) 



where (t) = '^p^p (''") creates an electron in the 
band state i:{r) = WpC'^''. The localizcid state is 
filled prior to turmeling, thus the density matrix in Eq.(5) 
is p = Peb+h. Here Pe = Ylpi^pa+ap + (1 - np)apa+) 
describes band electrons with energy distribution n(ep). 
The latter quantity can also be written as an exponential 
of a quadratic many-body operator, 



Pe = ^ exp ^- Apa+ap^ , 



e = 



1 - n(ep) 



with Z = Ylp{i + e'^p). One can bring G(r) to a 
standard form which depends only on the band electron 

variables by summing over the hole states [8,9]. This is 
achieved by disentangling b and 6+ from Eq.(5), 

ip+{-T)b{-T) = e-^^^Vi+Se*^^ = -Se-^'^i^^+e*^"^ 

and then using the commutation relations 

a+pe = e^ppea+ , a+e'^"^ = e-''^^ e^'^°^ a+ 

After summing over b, 6+ we obtain an expression 

G(t) = ^ M;,MpeVe-"p'^tr [e-'^'^e'^'>^ p^a+ap] (7) 
p,p' 

The central point of our approach is a relation be- 
tween the many-body operators in Eq.(7) and appropri- 
ate quantities (scattering operators and energy distribu- 
tion) defined in a single-particle Hilbert space. This re- 
lation holds [25] for any electron density matrix of the 
form of an exponential of a quadratic many-body opera- 
tor, such as Eq.(6). 



The advantage of introducing the single-particle scat- 
tering operators in the formalism at an early stage of 
the calculation is two-fold. First, wo bypass solution of 
the single particle scattering problem which requires re- 
summation of diagrammatic series [9] for band electron 
in the presence of a time-dependent scattering. Second, 
we shall be able to construct a non-perturbative solution 
applicable to an arbitrary energy distribution. 

We recall that Ref. [9] treats FES by solving the Dyson 
integral equation using a decomposition of quantities into 
analytic and anti-analytic functions of complex time vari- 
able, made possible by breaking the Hilbert space into the 
positive and negative frequency components. This ap- 
proach arises naturally in the equilibrium problem with 
a pure step n(e) but fails for a generic distribution. Below 
we develop an adequate replacement of this scheme. 

The discussion in the following two paragraphs closely 
follows that of Ref. [24]. First, we introduce an operator 
it) defined in the single particle Hilbert space of a band 
electron via the following operator product 



g-iHirgiWor^^ = ^"^ exp | ^ Wpyd+CLp, | (8) 

yP,P' 

The trace in Eq. (7) can be expressed through the oper- 
ator w as follows: 

tr {e-'^^^e''^°^pea+ap) = (i + 6-%^, det(i + e'^). (9) 

Our task is thus reduced to analyzing the quantity e*" 
which can be expressed through single particle operators: 



(10) 



with the single particle Hamiltonian operators /i^"'^) de- 
fined in Eq.(4) and {e~^)pp' = e~^''6pp'. These rela- 
tions help [25] to bring the determinant det (l -|- e*") to 

the form Z det ^1 — n(e) 

many-body effects fully accounted for by the algebra in- 
volved in the determinant construction. 

Next, the evolution operator product e^*'' '^e''* is 
related to the scattering matrix [24]. For one channel, 



^g-ih(i)Tgih(«)T^(g)j^ with the 



0<t,t' <T 

else 



(11) 



with the phase shift S = Si — Sq describing the effect of 
the resonance level changing occupancy. This gives 



det 



(1 + e*) = Z det (l + (^ - l)n) 



(12) 



Similarly, (l + e"'^) = {n{e) + (1 - n{e))S~^j n(e) 

where n and S are operators in the Hilbert space of func- 
tions of time. The Green's function (7) then becomes 
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G{t) = L{t)D{t), D = det (l + (5 - l)n) (13) 
L = ^u*u^e-'^'^(l-n(e)) (n + S-\l-h)y\ (14) 



with Ue = X^pWp'^(e — Cp)- The factors L and corre- 
spond, in the terminology of Rcf. [9] , to the open line and 
closed loop diagram contributions, respectively. 

Once the problem is reduced to analyzing certain one- 
particle operators there are two ways to proceed. Given 
that n is diagonal in the energy domain, while S is diag- 
onal in the time domain, one can choose cither represen- 
tation to analyze the quantities in Eq. (13). The former 
is convenient in equilibrium, since the T = Fermi dis- 
tribution is just a Cauchy kernel [9]. However, since for 

generic n(e) the kernel ht,t' = J n(e)e"^(*~* )(de/27r) is 
fairly complicated, the time representation does not ap- 
pear to be useful. Here, instead, we employ the energy 
representation. We note that the operator (5 — 1) has 
a double step structure 6{t)0{T — t) and argue that the 
contributions of the two steps can be treated as indepen- 
dent with logarithmic accuracy. For a single step in the 
time domain, the corresponding operator has the form of 
a Cauchy kernel in the energy dowMin. Such energy-time 
duality allows to perform calculation in essentially the 
same way as in the equilibrium problem, with the roles 
of energy and time interchanged. 

Since we are primarily interested in the power law 
exponent of G{t) rather than a prcfactor, let us con- 
sider replacing the double step 0{t)9{T — t) by a sum of 
almost nonoverlapping contributions 6{t)e~*/^ + 9{t — 
t)e~^'^~*^/'^ , t' < T. Such a replacement is reasonable 
since it preserves the steps a.t t = 0,t and thus af- 
fects the corresponding shakcup contributions merely by 
r changed to r' in the cutoff of the logarithms. (In addi- 
tion, we will have to adjust the extensive part InGun cx: r 
of the closed loop contribution as described below.) At 
the same time, since at t' < r the two terms do not over- 
lap, the operator quantities in Eq. (13) factor into two 
independent contributions. Employing this idea, we re- 
place the scattering operator S by a product T1T2, where 

{f,-l)t,t'=St,t'Xe{t)Ae-*^-', (15) 
(T2 - l)t,t' = St,t' X e{T - i)Ae(^-*)/"' (16) 

(A = e^*"^ - 1). We note that at t ~ 0, where fi time 
dependence has a step, the operator T2 is close to unity, 
while at t ~ T, where T2 has a step, Ti is close unity. 
This transformation allows to treat the contributions of 
Ti, T2 independently, which greatly facilitates analysis. 

We first analyze the open line contribution (14). In the 
S = T1T2 approximation, r' ^ r, the operator in Eq.(14) 
is factored into independent contributions as 



with Bj =t- — i. Let us write Bi in the energy domain: 



Bi = (e-2^^ - l)a, C7,,e 



27r e — e' — i/r 



^-7ZJ (18) 



Hence 1 + i?i(l - n) = I - a + af{e) with /(e) = 
j-g-2i5 _ _ n(e)) + 1. To invert this operator we use 
analytic properties of a. We note that &e,e' turns into a 
Cauchy kernel at large t'. Thus in this limit the opera- 
tor a projects to zero the functions F+ (e) analytic in the 
upper halfplane of complex e, Ime > 0, while a* = 1 — a 
projects to zero the functions Y-(e) analytic in the lower 
halfplane, Ime < 0. Conversely, a* operates as an iden- 
tity in the subspace of functions Y+{e), while a is an 
identity in the subspace of functions Y-{e). Hence it is 
convenient to employ factorization 

/(e)=r+(e)y-i(e) (19) 

The factors Y± are given in explicit form by 

lny±(e) = / ^^i#-de' (20) 

^ ' 2TriJe-€'±iO ^ ' 

Then, taking into account analytic properties of Y±, 
{l-a + afie))ll, = y_(e) ((1 - a)Yz\e') + aY^\e')) 

Similarly, the inverse -|- ^2(1 — n)^ is given by 

e-'^^F+ ^(e)((l - a*)Y+{e')+a*Y_{e'))e'''^ (21) 

where d*g, = (i/27r)(e — e'-|-i/r')~-^. After summing over 
e, e' in Eq.(14), we obtain L{t) = J2e L{e)e~^'^'^ where 



L{e) = \u,\\l-n{e))Y-\e) 



(22) 



To better understand this general result, let us consider 
the two-step distribution (1). Using Eq. (20), we obtain 



lnr+ = lln^-^ln^i^ 
TT /Ui — e TT /ii — e 



(23) 



Here jj.1^2 = Mi, 2 — i/T and 5 is defined by Eq.(3). Sub- 
stituting this result into Eq. (22), we obtain a split-peak 
structure with power law singularities at e = ^1,2: 



L(e) = 



|nep(l-n(e))(-eo)^^/- 
(e-/ii)2(5-^)A(e-M2)2^/" 



(24) 



At large |e| S> |/i2 — Mil; the power law form L oc e~'^^l'" 
matches the equilibrium result [9]. 

We now proceed to calculate the closed loop contribu- 
tion (13). First, consider a quasiclassical result, obtained 
by treating the time and energy as commuting variables: 



(n + 5-(i-n))"= n^(i + 4(i-n))"\ (17) In det (1 -n + 5n) = In (1 + MO) ^e' (25) 
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Thus we have det(l — n + Sn) = e~^'^ with complex 
C, = 7 + ieg, where the real part 7 describes broaden- 
ing of the FES singularity, while the imaginary part eo 
describes energy offset and can be absorbed in the phase 
factor e"*^"^. Evaluating the integral (25) for the two-step 
energy distribution (1), we obtain 

-7 = ^ In (1 - 4x(l - x) sin^ 5), 12 = H2 - Hi (26) 

Thus quasiclassical FES energy structure is a broadened 
step. The power law singularity appears only beyond the 
quasiclassical approximation. To describe it one has to 
account the contributions highly nonlocal in time, corre- 
sponding to many low energy particle-hole excitations. 

Again using the factorization approximation S = TiT2 
with soft cutoff e"*/"^ , we factor the determinant (14) as 



D = D1D2 = Yl det (1 + {fj - l)nj 



(27) 



j=l,2 



with the two factors accounting for the contributions of 
abrupt switching at t = and t = t. It is clear that the 
two determinants Di^2 are equal, therefore, it is sufficient 
to evaluate just one of them. Let us consider 

£>i = det(l + (fi - l)n) (28) 

The logarithm In D\ can be represented as a sum 



InDi = 2^iin + ^log 



(29) 



with Clin oc r and Ciog oc Inr. We have already es- 
timated the former (see Eq.(26)); to obtain the lat- 
ter, we consider variation AlnDi caused by a change 
in the distribution n{e). With the help of the formula 
Alndet;7 = tr (fZ-^Af/), we obtain 

Aln£)i = Ain \{l-u*+u*{Ah + l))~'^u*i^fi\ (30) 
To evaluate this expression we employ factorization 

An(e) + l=X_(e)X+i(e) (31) 
to rewrite the right hand side of Eq. (30) as 
tr (X+(e) [(1 - CT*)X+i(e') + a*XZ^{e')\ (T*A An(e')) 

Since -'^-i-(e) is analytic in the upper halfplane, /(I — 
(T*)e,e'X+^(e')<T*/ ^„de' = and Eq. (30) becomes 



Alnl?i =^tr {X+CT* XZ^ l^fi) 



(32) 



Prom Eq.(32) we extract the contribution corresponding 
to Clog: 



ACi, 



og 



I 

2^ 



/i 



ylAn(£) d 
+ An{e) de 



\nX+{e)de 



(33) 



So far the particular form of the distribution function 
n{e) did not matter, since Eqs.(22),(33) solve the prob- 
lem for arbitrary n(e). For the two-step distribution (1) 
the function lnX_|_(e) takes the following form 



\nX+{e) = In In — 

TT —C,0—l/T TV fli—t — l/T 



(34) 



with ^0 — Ep an ultraviolet cutoff. To find ACiog we 
consider variation of (1) with respect to x. 



Eq.(33) then yields a relation 



9Ci, 



og 



dx 



= iA 



2'k'^{1 + Ax) 



(35) 



Using the known value of Ciog for the equilibrium distri- 
bution, we integrate Eq. (35) over < a;' < a; and obtain 
the logarithmic term Ciog as a function of x, 



(1 + ^2^2)52/2^^ 



-(-ir^o)-''/'"' (36) 



Finally, restoring the exponential from Eqs. (25), (26), 
which is responsible for broadening of the FES, we obtain 



the closed loop factor D{t) 



exp(— 7t). 



The function D{e) ~ J eZ'^'^ D{T)dT defines broadening 
of singularities in the split FES, Eq.(2). We note that 
the broadening is relatively insignificant at J ^ 1, since 
ai,2 oc 6 while both 7 and the exponents in (36) are of 
order at small 5. Thus, although FES broadening is 
present for a split Fermi step, its magnitude is not large 
enough to smear the split peak FES profile. 
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